Multiplicative Levy noise in bistable systems 
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Stochastic motion in a bistable, periodically modulated potential is discussed. The system is 
stimulated by a white noise increments of which have a symmetric stable Levy distribution. The 
noise is multiplicative: its intensity depends on the process variable like |a;|~*. The Stratonovich 
and Ito interpretations of the stochastic integral are taken into account. The mean first passage time 
is calculated as a function of 6 for different values of the stability index a and size of the barrier. 
Dependence of the output amplitude on the noise intensity reveals a pattern typical for the stochastic 
resonance. Properties of the resonance as a function of a, 6 and size of the barrier are discussed. 
Both height and position of the peak strongly depends on 6 and on a specific interpretation of the 
\ stochastic integral. 

I. INTRODUCTION 

Ubiquity of the normal distribution in nature is a consequence of its stability. However, stable distributions are not 
restricted to the Gaussians: processes which involve discontinuous trajectories and algebraic tails ^ where 
O, 0<a<2isa stability index, can be stable as well. Divergent moments of those distributions reflect a presence 
^ • of long jumps (Levy flights). The normal distribution corresponds to a = 2. The Levy stable statistics with a < 2 
is important when one considers realistic phenomena, where complexity, non-uniformity and long-range correlations 
play a role. For that reason, they are encountered and studied, for example, in biology sociology [1| and finance 
^ \ [3, |j| . A characteristic feature of many problems in those fields is a non- homogeneous structure of the medium which 
c/2 ■ results in an anomalous transport, as it is the case for porous, disordered materials [5^^ and folded polymers [6j. The 
non-homogeneity of the medium has significant consequences for a stochastic description of such systems in terms of 
the Langevin equation; the noise can no longer be additive and must include a dependence on the process variable. 
Despite that, multiplicative problems with Levy flights are rarely discussed in the literature. 

It is a well-known fact that the Langevin equation with the multiplicative white noise is not unique and requires an 
(— ( ■ additional interpretation of the stochastic integral. Two of them are of particular interest: the Ito interpretation (II), 
Q '[ for which the integrand is evaluated before the noise action, and the Stratonovich one (SI) which takes into account 
O ■ the process value both before and after the stochastic stimulation. SI is distinguished since it constitutes a limit 
of the coloured noise in a form of the fractional Ornstein-Uhlenbeck process [7]. On the other hand, II applies for 
, discrete systems and systems with a large inertia [tI, HI • II leads to a fractional Fokker-Planck equation with a variable 
^ ■ diffusion coefficient [§| and it differs from that for SI only by a spurious drift if a = 2. For non-Gaussian processes a 
, relation between both interpretations is more complicated and predictions qualitatively different: for SI, in contrast 

■ to II, slope of the distribution tail is modified by the multiplicative factor [l^, Ull- In particular - in the case of the 
' linear driving force - that factor may sufficiently weaken the noise to make the variance of the output signal finite. 
. Slopes corresponding to the finite inertia lie between values for II and SI 0] which emphasises importance of taking 

CN| ' into account both interpretations. 

, One also can achieve the finite variance by introducing a nonlinear deterministic force to a system with the additive 
noise From the physical point of view, the double- well potential is of particular interest since it allows us to 

" . [ , study confined systems and, in particular, to model the barrier penetration. The double-well system was examined 
^ ■ predominantly for the Gaussian case [l3l ; then the first passage time problem resolves itself to a dynamical equation 
with ordinary boundary conditions [14{. The barrier penetration for the Levy flights was considered in Ref . [l^ [l^ . 
New effects emerge when one introduces - in addition to the double-well potential - a time-dependent periodic force. 

■ A rate of the jumping between the wells depends nonlinearly on the noise intensity and then a tuning of that intensity 
yields a matching of two characteristic frequencies: the noise-induced rate and a frequency of the deterministic, 
periodic force. That phenomenon, called a stochastic resonance [13], is well-known for the Gaussian noise [l^ and 
was analysed also for the multiplicative case [l^ . If the additive noise involves the Levy flights the stochastic resonance 
also emerges; it was observed in systems with both the double- [l^, [2l| and single- well potential [22j. 

In this paper, we address a problem of the resonant switching between potential wells for the case of the multiplicative 
Levy-stable stochastic stimulation. We demonstrate to what extend properties of the system are modified, in respect 
to the known cases, when one introduces long jumps and takes into account a dependence of the noise on position 
for both SI and II. In Sec. II the system, which includes the bistable potential and the oscillatory force, is defined and 
a mean first passage time calculated. Existence of the stochastic resonance is demonstrated, and dependence of the 
signal amplitude on the system parameters discussed, in Sec. III. 
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II. PROBABILITY DISTRIBUTIONS AND MEAN FIRST PASSAGE TIME 



A general stable Levy distribution is expressed in terms of three parameters: the stability index a (0 < a < 2), 
which determines the tail of the distribution, the asymmetry parameter (skewness) and the translation parameter. 
Since, in the following, we restrict our analysis to the symmetric distributions, the index a completely determines the 
driving noise. The characteristic function for this case reads 

p^{k) = exp{-a"\kn, (1) 

where a measures the apparent width of the distribution. If a = 2 Eq.(IT]) represents the normal distribution, otherwise 
the tails are algebraic, ~ |^|^^~", and the stochastic trajectory is discontinuous in a sense of the Lindeberg condition 
[23 |. The variance is infinite, as well as all higher moments. 
We consider a motion in the time-dependent bistable potential 

V{x,t) ^ -^x'^ + ^x* ~ Aoxcos{ujot) (2) 

where the amplitude Aq and the frequency ujq are constant. The system is driven by the noise G{x)rj{t) and we 
assume that increments of the noise rjit) are distributed according to the characteristic function ([T|). Then the 
Langevin equation is of the form 

x^ -dV{x,t)/dx + G{x)f]{t). (3) 

Intensity of the noise algebraically depends on the position, 

G{x) = \xr\ (4) 

where 9 > —1. In this paper we take into account two interpretations of the multiplicative noise in Eq.Q: the 
Ito interpretation and the Stratonovich interpretation. Components of the discretized stochastic integral in II are 
defined by G[x{ti-i)\[T]{ti) — 77(^^-1)], i.e. G[x{t)] is calculated before the noise acts. In SI, in turn, one takes a middle 
point: G[{x{ti-{) + x{ti)) /2][r]{ti) — 77(^^-1)]. A difference between both interpretations, as well as probability density 
distributions for some nonlinear potentials, is well-known for the Gaussian processes [13, HI]. The Fokker-Planck 
equation, corresponding to Eq.Q in II, involves a fractional derivative and a variable diffusion coefficient ,9]: 

^p(a:, i) = - ^n^)Pi^, t) + a" ^ |G(x) t) , (5) 

where F{x) — —dV{x,t)/dx. 

A technical advantage of SI consists in an admissibility of the variable change which allows us, in the one-dimensional 
case, to reduce Eq.([3]) to the Langevin equation with the additive noise. Validity of the ordinary rules of the calculus 



was rigorously proved for a = 2 23|. It was confirmed for a < 2 by the numerical simulations but the linear case 
required introduction of a cut-off to the distribution [l^, UH. On the other hand, rules of the ordinary calculus 
apply to systems with the coloured noise and the white-noise limit exists which was demonstrated for the generalised 
Ornstein-Uhlenbeck process . After the transformation 



the Langevin equation takes the form 



where the effective potential reads 



y^-dV{y,t)/dy + f^it) (7) 



Wells of the time-independent part of the potential are positioned at ym = ±(a/6)*^^+^''/^/(l + 6) and their depth 
declines with 6. 

Properties of the stochastic system are different for both interpretations of Eq.Q. The case V{x,t) = was 
discussed in Ref . [l0| . For II, the distribution tail does not depend on - it falls like |a;|~^~" - whereas SI leads to the 
dependence \x\~^~°'~^°' . As a consequence, variance can be finite for SI. Result for the potential V{x,t) = —Xx'^/2, 
where A > 0, is similar. 
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FIG. 1: (Colour online) The probability density distributions for II (blue dotted lines with symbols) and SI (red solid lines with 
symbols) calculated with a — 1.5 at t = 2. The case 9 = 2 is marked by squares and 9 — —0.5 by circles. Other parameters 
are a = 4, b = 1, Aq = 3, uq = 5 and a = 10. The case ^ = with a = 1.5 and 2 are marked by the solid and dotted lines, 
respectively. 



Presence of the nonlinear force modifies tails of the distribution, compared to linear case, making them steeper 
because, despite long jumps, the potential well prevents the particle from escaping to large distances. As a consequence, 
the variance may become finite. It is always the case for a quartic oscillator both with the additive noise [13] and the 
multiplicative noise in SI if > —1 [ll|. Motion in the double- well potential with the a— stable additive driving was 
considered by Ditlevsen [isj . The waiting and barrier penetration time - defined in terms of the reflection and/or 
absorbing barriers - is of particular interest for that system. An influence of the algebraic tails of the distribution for 
a < 2 on the average transition time between the wells becomes visible if width of the barrier is sufficiently large; then 
the transition may be caused by a single jump [l5j|. Otherwise the motion is dominated by a continuous superposition 
of many small jumps, similarly to the Gaussian case, and the transition time depends exponentially on the potential 
depth according to the Arrhenius formula. The mean first passage time displays a bell-shape as a function of a: it 
rises with a up to a = 1 and then falls [l6j . In the case of the multiplicative noise, a behaviour of the transition time 
as a function of system parameters depends on the specific interpretation of the stochastic integral . 



A further generalisation of the above system includes the periodic driving force. For this system, Eq.([3]), we define 
the mean first passage time Tp as the average time needed to pass for the first time from the right minimum of the 
potential to the left one. Therefore, we introduce an absorbing barrier at a; = Xm — —\/oJb in a sense that trajectories 
for which x{t) < Xm are absorbed. The fractional Fokker-Planck equation corresponding to Eq.(|31) involves a nonlocal 
boundary conditioUjSince the jumping particle may skip the boundary without hitting it, and the problem cannot be 
solved analytically |lQ]. In the present paper, the stochastic trajectories were calculated by the numerical integration 
of Eq.([31). In the case of II, a simple Newtonian method was applied. For SI, in turn, we first transformed the variable 
X ^ y, according to Eq.®, and then applied the Heun method to the equation with the additive noise, Eq©: 

Vi+i = yi + hF{y„ti) (9) 
y^+l = y^ + ^[F{y^,U)+Fiy,+uU+l)]+h^/''^, 



where F{y,t) — -~dV{y,t)/dy. 
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FIG. 2: (Colour online) Mean passage time from the left to the right well of the potential for II (black dotted lines) and SI 
(red solid lines) as a function of 6. The case a = 2 is marked by squares and a — 1.5 by circles. Other parameters are a = 4 
(full symbols), a — 1 (empty symbols), b = 1, Ao = 3, uJo = 5 and a — 10. 



The multiplicative term G{x) qualitatively modifies the distribution. First, we consider the case without the 
absorbing barrier. The distribution possesses two peaks, at negative and positive x, corresponding to the well positions. 
Their height decreases with a and, for 6 < 0, they vanish altogether. Instead, a peak at the origin emerges as a result 
of increased contribution from the interwell motion. Fig.l presents the probability density distributions obtained by 
a numerical integration of Eq.Q with the initial condition p{x, 0) = S{x — Xm)- A relatively large noise intensity has 
been appHed and the interwell motion prevails. In contrast to the case 9 — 0, the distribution is strongly anisotropic 
for both interpretations. Since the noise vanishes at the origin for 9 < 0, the particle abides near the top of the 
barrier. On the other hand, for the case 9 > both peaks are present even for a large a. Those conclusions apply for 
both II and SI. The case 9 = (additive noise) exhibits small maxima at the well positions for a = 1.5 whereas the 
Gaussian case is almost uniform in a wide range of the argument. 

The mean passage time as a function of 9 for various configurations is presented in Fig. 2. Results depend both 
on the stability index a and on the particular interpretation of the stochastic integral. Moreover, a sensitivity on 
the potential geometry is visible. First, let us consider the case of the relatively wide barrier (a = 4). Depth of 
the effective potential for SI, Eq.([5]), diminishes with 9 and one could expect a smaller Tp. However, we know from 
Ref.[lB| that, for a system with the additive noise, Tp depends rather on the width of the barrier than on the potential 
depth if jumps determine the barrier penetration. We observe a similar effect: position of the effective well reaches 
its minimum at 9m = 2/ln(a/6) — 1 which value corresponds to the minimum of Tp in the figure. Nevertheless, there 
is no qualitative difference between the cases a = 1.5 and 2. Shape of the curves Tp{9) for II is similar but Tp rises 
faster for > 0, as a result of a strong noise damping. The passage time as a function of 9 for a = 1 is different than 
for a = 4 because then ym{9) = 1/(1 + 6), 9m does not exist and Tp{9) monotonically falls. Since 9 > —1, such a 
behaviour is expected for all narrow barriers, namely if a < 6. Obviously, Tp in Fig. 2 decreases with a because the 
noise intensity cr". 



The rate of jumping between potential wells due to the noise can match the frequency wq of the deterministic 
oscillatory force revealing a resonant pattern. This phenomenon has been extensively studied for the Gaussian 



III. STOCHASTIC RESONANCE 
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FIG. 3: (Colour online) The probability density distributions for both interpretations calculated ai t = 2 with the parameters: 
a = 1.5, a = 4, 6 = 1, ^0 = 3, ojo = 5 and a = 2. Cases with the multiplicative noise were calculated with 9 = —0.5. The red 
solid lines are of the form x~'^ with fi = 5.5 and 4.75. 

processes and observed also for the general Levy statistics of the noise. It was demonstrated in Ref.flOl that the 
signal-to-noise ratio exhibits a resonant amplification for many systems characterised by the long jumps. The spectral 
amplification and dynamic hysteresis loops were calculated for the asymmetric Levy noise with various a and skewness 
parameter ^21]; the amplification dwindles with decreasing a. On the other hand, the stochastic resonance phenomenon 
for a linear multiplicative Gaussian noise was studied in Ref. (19j and attributed to the transition from an intrawell to 
an interwell modulated dynamics. 

A useful quantity to describe a stochastic resonance is a correlation between the output process x{t) and the periodic 
driving, x. It equals the first Fourier coefficient which follows from the expression [26| 



where T — 2tt /loq and p{x, t) results from Eq.Q. The density t) is asymptotically periodic and an averaging over 
the period produces a stationary distribution. The spectral amplification at the frequency cjq is given by a ratio of 
the integrated power stored in spikes of the power spectrum to the total power carried by the oscillatory force, 77 = 
4(x/Ao)^ where X = ReMi [11]. The signal-to-noise ratio (SNR) also can be expressed by Mi : SNR= Ait\Mi\^ / S%{ujo) 
where S'^^ is the total power spectral density. The latter quantity is directly related to the autocorrelation function 
and its existence requires the finite variance. That condition is satisfied for systems with the quartic potential in SI 
since then the tail ^ \x\^^^°'^°'^ On the other hand, an analysis for both the free particle and the linear force 
indicates that in II shape of the distribution tails does not depend on 9. Fig. 3 shows that also for the general system 
^ the tails fall faster than for the simple Levy motion: the slope rises with 9 for SI whereas it does not depend on 
9 for II. Therefore variance remains finite when Levy flights are included and the spectral quantities, r] and SNR, are 
finite and well determined. 

For small amplitudes and large time, a response of the system to the periodic stimulation is proportional to x, 
{x{t))as = xcos{uJot — (f>), where is a phase lag [l^. The stochastic resonance is usually detected by plotting x as 
a function of the stochastic stimulation a. That function does not monotonically diminish but exhibits a bell-shaped 
structure maximum of which corresponds to a matching condition of two frequencies: uja and a rate of the stochastic 
jumping, 1/Tp. Then we can expect that a position of the resonance is directly related to the mean passage time and 
satisfies the approximate condition 




(10) 



uJoTp{a,a,9) 1. 



(11) 
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FIG. 4: (Colour online) a; as a function of cr for 9 — 1.1 (solid lines) and —0.6 (dotted lines). The following values of a 
are presented: 2 (black squares), 1.8 (blue circles), 1.5 (green triangles) and 1.2 (red diamonds). Lines with empty symbols 
corresponds to 6 — 0. Other parameters are a = 4, 6 = 1, = 3, uo = 5. Averaging was performed over 2 x 10^ events. 
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FIG. 5: (Colour online) Left part: ar versus 9 for SI (solid lines) and II (dotted lines). The case o = 4 was calculated with 
a = 2 (black squares) and 1.5 (green circles), the case a = 1 with a = 1.8 (blue stars). Other parameters are the same as in 
Fig. 4. Right part: height of the peak versus 9. 
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FIG. 6: (Colour online) The same as Fig. 4 but for a = 1 and a — 2, 1.8 and 1.5. 

Hight of the peak rises with the amplitude Aq and diminishes with the frequency ujq. In the following, we discuss an 
emergence of the stochastic resonance in the system ([3]) and demonstrate how its position and intensity depends on 
a, 0, as well as geometry of the potential. 



Fig. 4 presents results of the calculations: a dependence of x on tr for two values of 9, both positive and negative; 
comparison with the additive noise is included. The resonance structure is visible for all the presented cases. Re- 
markable differences emerge when one compares both interpretations. For SI, x is large if < and a position of the 
peak, ar, occurs at larger noise intensity than for > since then the effective barrier is wider. The resonance for 
the positive 6 emerges at roughly the same a as for 9 = 0. The peak is much lower for II if < and corresponds to 
a smaller a than for > 0. The latter case, governed by Eq.®, is characterised by a weak noise intensity when the 
particle remains far from the origin and that effect must be compensated by a larger a. Dependence on a is similar 
for both interpretations: ar shifts towards smaller a and rises with a for 9 < 0. Ur falls with a also for 

9 > Q but this effect is much weaker, especially for II. A dependence of (Jr on 9 is presented in Fig. 5 for a = 1.5 
and 2. The difference between both interpretations is especially pronounced for 9 < Q when Ur remains constant for 
II, in contrast to SI. For > 0, all presented cases exhibit a monotonically increasing function (Jr{9)- Height of the 
resonance, x^, is also presented in Fig. 5 as a function of 9. This dependence is very weak for II, in contrast to SI 
for which we observe a pattern similar to Urif)) with a minimum at the same point. Results for both values of a are 
almost the same. An apparent similarity between Fig. 2 and 5, which is especially striking for SI, is a consequence of 
the condition (jlip : since the curves Tp{9) for different <j are, approximately, parallel to each other, they mimic the 
dependence Ur{9). 

A change of the potential shape qualitatively modifies Tp and one can expect differences also for the stochastic 
resonance. This conclusion is illustrated in Fig. 5 which compares the case of the wide barrier (a = 4) with the case 
of the narrow barrier (a = 1) for both interpretations. Though results for II do not differ substantially - Ur rises for 
both geometries - the curve for SI and a = 1 is flat resembling a similar behaviour of Tp, cf. Fig. 2. Comparison of the 
peak heights is also presented in Fig. 5. We conclude that the amplitude is not sensitive not only to a but also to the 
barrier width, in contrast to peaks' position. The amplification is much stronger for SI than for II and the dependence 
Xr{9) is determined mainly by a specific interpretation of the stochastic integral. Curves x{a) for the narrow barrier, 
shown in Fig. 6, are similar to the case a = 4 (Fig. 4); difference resolves itself to a smaller cr^ for SI if > 0. 
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IV. CONCLUSIONS 

Properties of systems containing stochastic force, which is characterised by the stable distribution with the heavy 
tails, differ from those for the normal distribution possessing, in particular, the infinite variance. Systems with the 
multiplicative noise exhibit a more complicated behaviour and the variance may be finite, despite the algebraic tails. 
Nonlinear potentials, in turn, makes the distribution tails steeper and may invoke an effective trapping. In such 
systems motion is limited in space and asymptotics of the driving noise plays a relatively minor role, compared to the 
diffusion problems. Motion in the double- well potential is a simple example. 

We have discussed such a system: the anharmonic oscillator with the periodic modulation and the stochastic 
stimulation in a form of the stable, multiplicative process. The multiplicative factor, assumed in the algebraic form, 
makes the distribution anisotropic; the point x — becomes an attractor for 6 < since noise vanishes at the origin. 
The mean first passage time depends on a particular interpretation of the stochastic integral; for SI it is determined 
rather by the size of the effective potential than by the effective well depth. That conclusion, however, is valid only for 
a relatively wide barrier and then the geometry is important for the dependence Tp{9). The passage time smoothly 
rises with a for both geometries and no qualitative difference in respect to a = 2 has been observed. 

The system ([3]) possesses two characteristic frequencies, Wq and 1/Tp, and, when they match, the resonance occurs. 
The resonance shape of the dependence xia) has been observed for all analysed cases. Since the variance is finite 
also for processes with the Levy flights, both the spectral amplification 77 and SNR are well determined. For SI, 
height of the resonance is sensitive on 9 and the amplification rapidly rises when becomes negative. Position of 
the resonance, in turn, mimics the dependence Tp(d). Therefore, size of the barrier is important for the resonance 
position: ct^ is almost constant for SI if a = 1, in contrast to the case a = 4. We finally conclude that - though height 
and position of the peaks depend on a - the bistable system with the Levy flights has similar properties as the case 
of the Gaussian processes. On the other hand, dependence of the noise on the process value essentially modifies - 
differently for each interpretation of the stochastic integral - the distributions, the passage time and properties of the 
stochastic resonance. 
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